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O ’ Abstract 

A sensitivity  study  was  made  on  four  variables  which 
must  be  specified  to  ccmipute  the  thermal  radiation  fluence 
in  the  atmosphere  from  a nuclear  burst.  These  parameters 
are;  the  choice  and  number  of  energy  bands  used  to  specify 
the  atmospheric  attenuation  coefficients,  the  relative  lm> 
portance  of  scatter,  the  significance  of  source  (fireball) 
temperature  and  the  Importance  of  height  of  burst.  These 
parameters  were  examined  by  computing  the  thermal  fluence 
as  the  variables  were  changed,  one  by  one.  The  number  of 
groups  was  varied  fr(xn  four  to  sixteen.  The  source  tempera, 
ture  was  examined  over  a range  from  3500^  to  8500®K  and 
(j  burst  altitudes  to  30  km  were  considered.  The  Importance 

of  scatter  was  evaluated  by  computing  the  fluence  first 
with  an  attenuation  factor  composed  only  of  absorptl(»i  and 
then  with  an  attenuation  which  was  the  sum  of  absorption 
and  scatter.  All  of  these  cal^'ulatlons  were  carried  out 
with  a POBTBAN  IV  computer  code  prepared  by  the  author.  The 
code  Is  described  here  for  possible  future  use.  The  results 
of  the  parametric  variation  Indicated  that:  one  needs  14 

energy  bands  In  the  Infrared  region,  scatter  Is  Important  In 
computing  fluence  at  altitudes  below  five  km,  the  computed 
fluence  Is  relatively  Insensitive  to  source  temperature 
over  the  ra^e  considered  and  that  burst  height  changes  can 
be  modeled  by  changing  source  temperature. 
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A SENSITIVITY  STUDY  OP 
THEBHAL  BADIATION  PLUENCE 
PROM  A NUCLEAR  AIR  BURST 

I.  Introduction 

Calculations  of  the  effects  of  nuclear  weapons  are 
Important  In  predicting  the  survivability/vulnerahility  of 
Air  Perce  systems.  Thermal  radiation,  one  of  the  major 
effects  generated  by  the  fireball  of  a nuclear  explosion, 
contributes  approximately  30  to  40  i>ercent  of  the  total 
weapon  yield  from  an  air  burst  (Ref  2:7)# 

The  objective  of  this  thesis  is  to  present  a sensitivity 
study  of  the  thermal  fluence  from  a nuclear  air  burst  in- 
cident on  a receiver  in  the  atmosphere,  A nuclear  air 
burst  is  defined  as  a burst  at  an  altitude  of  less  than 
100,000  ft,  or  about  30  km.  Por  the  discussion  and  examples 
to  follovr,  both  burst  and  receiver  are  equal  to  or  less 
than  30  km.  Virgin  radiation  is  that  fraction  of  the  thermal 
radiation  which  is  unscattered  and  not  absorbed  by  the  at- 
mosphere. When  calculating  the  virgin  fluence,  all  scatter 
emd  absorption  results  in  a loss  of  radiation  to  the  re- 
ceiver, Pluence  is  the  time  Integrated  energy  flux  and  is 
measured  in  units  of  energy  per  unit  area. 
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Chapter  II  of  this  report  presents  a short  overview 
of  the  theory  of  thermal  radiation  frcmi  an  air  burst. 

Here,  the  source  of  the  radiation  Is  presented  along  with 
a transmission  model  based  on  diffusion  theory,  In  order 
to  develop  a physical  basis  for  the  calculation  of  fluence. 
The  energy  range  of  the  thermal  photons  and  the  possible 
mechanisms  of  photon  Interaction  with  the  atmosphere  are 
examined.  The  degree  of  photon  Interaction  with  the  atmos- 
phere Is  measured  by  the  photon  cross  sections  determined 
In  Chapter  III.  Atmospheric  models  of  the  earth's  atmos- 
phere are  compared  and  cross  sections  are  determined  by 
the  use  of  transmittance  data  developed  by  B.  A.  HcClatchey 
of  the  Air  Force  Cambridge  Besearch  Laboratory,  Chapter  IV 
presents  a series  of  sensitivity  studies  made  by  calculating 
thermal  fluence  from  an  air  burst.  The  parameters  studied 
are:  the  choice  and  number  of  energy  bands  used  to  specify 

the  atmospheric  attenuation  coefficients,  the  relative  Im- 
portance of  scatter,  the  slgniricance  of  source  (fireball) 
temperature,  and  the  Importance  of  height  of  burst. 

All  the  cross  section  determinations  and  sensitivity 
studies  are  combined  In  the  development  of  a thermal  code 
In  Chapter  V.  The  results  of  computer  problems  with  200- 
klloton  bursts  summarize  the  conclusions  of  the  sensitivity - 
study.  The  code  Is  developed  using  basic  POBTBAN  IV  pro- 
gramming language  which  can  be  followed  by  a student  with 
cnly  a bssie  'Understanding  of  POHTEAN,-- 


one/ph/75-7 


I V.  ’ II,  Theory  of  Thermal  Radiation  Phenomenon 

i 

! 

I Source 

The  source  of  thermal  radiation  In  a nuclear  explosion 
Is  the  fireball.  In  an  air  burst,  the  emission  of  thermal 
radiation  Is  In  two  pulses  which  correspond  to  the  rise 
and  fall  of  the  apparent  surface  temperacure  of  the  fire- 
ball. The  two  pulses  are  Illustrated  In  Figure  1. 


TIME  AFTER  EXPLOSION 
(RELATIVE  SCALE) 


Fig.  1.  Emission  of  Thermal  Radiation  In  Two  Pulses  In 
an  Air  Burst  (Ref  6:45). 

In  the  first  pulse,  there  Is  only  about  one  percent  of  the 
thermal  radiation  because  of  Its  short  duration,  which  Is 
about  0.1  second  for  a one-megaton  explosion.  The  second 
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pulse  may  last  for  several  seconds  and  carries  about  99 
percent  of  the  total  thermal  raalatlon  (Ref  6:45) • The 
first  pulse  is  neglected  In  calculations  of  thermal  radla. 
tlon  from  an  air  burst. 

The  average  surface  temperature  of  the  fireball  for  the 
second  pulse  Is  approximately  6000^K.  The  fireball  acts 
much  like  the  sun  and  can  be  approximated  by  a blackbody 
radiator,  so  a 6000^K  blackbody  Is  often  used  as  the  source 
In  a simple  model  of  a nuclear  air  burst.  The  nuclear  air 
burst  spectrum  will  vary  from  the  6000^K  blackbody  mostly 
In  the  ultraviolet  region  (Ref  9). 

Transmission 

The  theory  of  transmission  of  thermal  radiation  can 
be  developed  from  the  diffusion  equation.  Assume  a point 
source  with  a homogeneous  atmosphere  and  all  scatter  as 
removal.  Because  of  the  absence  of  scatter,  the  thermal 
radiation  acts  like  streaming  particles.  Therefore  the 
fluence,  F(r),  Is  equal  to  the  net  current  and  Is  a function 
of  the  distance  from  the  source.  Under  these  conditions, 
the  dlffuplon  equation  can  be  written 

3^  + ip(p)  + iiP{r)  - 0 ,,, 

dr  r (1) 

where  r Is  the  distance  from  the  source  and  (x  Is  the  total 
cross  section  (per  unit  length),  which  is  defined  as  bhe 
absorption  plus  scatter  cross  sections.  Rearranging 
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equatlcm  (1)  and  Integrating  from  the  radius  of  the  fireball, 
^fb*  receiver,  B,  gives 


B 


dP(r) 


/r 

-'Tfb 


/2dr./, 


dr 

fb 


which  Is 


rLP(r^^^) 

P(H)  « exp  (-u(B  - r^^,)) 


Applying  boundary  conditions 


lim  (4irr%(r))  » f 


'-►'■fb 

where  I Is  the  thermal  yield  (energy) 

P(r,b>  - ^ 


4tr; 


fb 


(2) 


(3) 


(^) 


(5) 


and  substituting  Equation  (5)  Into  Equation  (3)  yields 


P(B)  - -I—  exp(-pB)  .. 

4vB'^  (o) 

V 

where  the  fireball  Is  assumed  the  point  source  such  that 
'•fb  Is  zero.  The  expression  exp.(-pB)  Is  the  atmospheric 
transmittance,  T.  Atmospheric  transmittance,  through  the 
cross  section  term  p,  is  a function  of  many  variables;  wave- 
length of  radiation,  path  length  of  transmission,  atmospheric 
gases,  pressure,  temperature,  amounts  of  aerosols  and  the 
size  distribution  of  the  aerosols  (Ref  4:81).  If  a point 
source  Is  assumed  cmd  effects  of  the  atmosphere  are  Ignored, 
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the  thermal  fluence  at  the  receiver  could  be  calculated  by 


Y 


(7) 


where  P Is  the  thermal  fluence  (energy  per  unit  area),  and 
Y Is  the  thermal  yield  (energy),  and  B is  the  path  length 
between  burst  and  receiver  (length).  The  only  attenuation 
of  thermal  radiation  would  be  spherical  divergence  given 
by  47rH^. 


Spectrum  from  an  Air  Burst 

The  range  of  thermal  photon  energies  from  the  source 
Is  assumed  to  correspond  to  a 6000°K  Plancklan  distribution. 
Using  BAND  Plancklan  tables  (Ref  5)»  the  range  from  0.2  to 
4.0  microns  contained  approximately  99  percent  of  the 
thermal  radiation  emitted.  The  thermal  spectrum  of  0.2  to 
4.0  microns  can  be  divided  Into  three  broad  bands;  ultra, 
violet,  visible  and  Infrared^  Figure  2 gives  the  band 
limits  In  terms  of  wavelength,  frequency  and  energy. 

In  terms  of  percent  of  thermal  yield  In  each  band, 
approximately  13.8  percent  Is  In  the  ultraviolet,  37.5  per- 
cent In  the  visible  and  47.6  percent  In  the  Infrared. 


Photon  Interactions  with  the  Atmosphere 

The  total  cross  section,  |i,  determines  the  degree  of 
Interaction  between  the  thermal  radiation  and  the  atmosphere. 
The  atmosphere  consists  of  molecules  and  aerosols,  hole- 
cules  are  very  small,  on  the  order  of  10"^  to  10"®  microns. 
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two  mechemlsms  of  Interaction;  absorption  and  scatter* 

Absorption.  The  absorption  cross  sections  for  molecules 
are  a function  of  not  only  the  amount  of  absorbing  gas,  but 
also  the  local  temperature  and  pressure  of  the  gas.  The 
variation  of  the  molecular  absorption  with  wavelength  I3 
very  complicated,  being  a highly  variable  function  of 
wavelength  ( Ref  8:11). 

The  absorption  cross  sections  for  aerosols  depend  on 
the  number  density  and  size  distribution  of  aerosols  as 
well  as  the  complex  Index  of  refraction.  The  amount  of 
absorption  can  be  derived  theoretically  for  spherical 
particles  by  use  of  Mle  theory.  The  absorption  cross 
section  due  to  aerosol  distributions  Is  a smooth  function 
of  wavelength. 

Scatter.  The  scatter  cross  sections  for  molecules 
depend  only  on  the  number  density  of  molecules  In  the 
radiation  path.  The  wavelength  dependence  of  molecular 
scattering  Is  very  nearly  proportional  to  tne  Inverse 
fourth  power  of  wavelength.  All  atmospheric  molecules 
scatter  radiation.  If  the  atmospheric  molecules  are  assumed 
to  be  spherical  and  small  relative  to  the  wavelength  of  the 
Incident  radiation,  the  phenomenon  Is  called  Rayleigh 
scattering. 
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The  scatter  cross  section  for  aerosols  depends  on  the 
number  density  and  size  distribution  of  aerosols.  The 
amount  of  scatter  can  be  derived  by  use  of  Hie  theory  and 
is  a smooth  function  of  wavelength. 

The  determination  of  atmospheric  attenuation  is  reduced 
to  a calculation  of  the  attenuation  due  to  aerosol  and 
molecular  cross  sections.  The  total  aerosol  cross  section 
is  the  sum  of  the  scatter  and  absorption: 

"“a  ■ «.  + fa  (8) 

where  a.  is  the  aerosol  absorption  cross  section  (per  unit 

cL 

length)  and  is  the  aerosol  scattering  cross  section  (per 
unit  length) . Likewise,  the  molecular  cross  section  is  the 
sum  of  the  molecular  scattering  and  absorption: 


a 


m 


a_  + 
m 


(9) 


where  is  the  molecular  absorption  cross  section  (per  unit 
length)  and  is  the  molecular  scattering  cross  section 
(per  unit  length). 

The  total  cross  section,  is  the  sum  of  the  molecular 
and  aerosol  cross  sections  and  is  often  called  an  attenuation 
cross  section. 
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III.  DeteiTnlnation  of  Photon  Cross  Sections 


Atmospheric  Models 

The  atmospheric  models  considered  to  determine  cross 
sections  used  In  the  calculation  of  virgin  thermal  fluence 
were  models  compiled  by  Air  Force  Cambridge  laboratory 
(APCBL)(Ref  8).  Six  model  atmospheres  are  defined  for 
temperature,  pressure  and  absorbing  gas  concentrations.  In 
addition  to  these  atmospheres,  two  aerosol  models  have  been 
defined  and  are  available  for  use  with  the  six  models. 

Five  of  the  model  atmospheres  compiled  by  AFCBL  repre- 
sent different  climatological  conditions  that  are  typical 
for  seasonal  conditions  and  geographical  locations,  namely. 
Tropical,  Midlatitude  Summer,  Midlatitude  Winter,  Subarctic 
Summer  and  Subarctic  Winter.  The  sixth,  the  U.  S.  Standard 
Atmosphere,  1962,  was  also  compiled,  with  an  appropriate 
water  vapor  and  ozone  distribution.  The  U.  S.  Standard 
Atmosphere  provides  a single  appropriate  atmospheric  model 
for  comparison  within  the  Atmospheric  Physics  community. 

The  two  aerosol  models  compiled  by  AFCBL  describe  a 
"clear"  and  "hazy"  atmosphere  corresponding  to  a visibility 
of  23  And  5 hm  respectively  at  ground  level.  The  aerosol 
size  distribution  function  for  both  models  Is  the  same  at 
all  altitudes.  The  aerosols  considered  vary  In  size  from 
0.02  to  10  microns.  The  "hazy"  model  Is  Identical  to  the 
"clear"  model  for  altitudes  equal  to  or  greater  than  five  km. 
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For  altitudes  less  than  five  km,  the  total  number  of 
aerosol  particles  Increases  exponentially  to  reach  a value 
at  ground  level  corresponding  to  a ground  visibility  of 
approximately  five  km. 

For  the  purpose  of  this  study,  the  U.  S.  Standard 
Atmosphere  and  the  Tropical  Atmospheres  coupled  with  the 
"clear*  and  "hazy"  aerosol  distributions  were  used  to  deter- 
mine cross  sections  for  the  calculation  of  virgin  fluence. 

The  U.  S.  Standard  Atmosphere  was  chosen  to  represent  a 
"dry*  atmosphere  and  the  Tropical  Atmosphere  a "wet*  atmos- 
phere. Tables  I and  II  give  the  concentrations  and  meteoro- 
logical condition  as  they  vary  with  height  for  the  atmospheres 
chosen. 

Transmission  Data 

Cross  sections  were  determined  by  the  use  of  transmission 
data  compiled  by  McClatchey  of  APCBL  (Bef  8).  The  U.  S.  Stan- 
dard and  Tropical  Atmospheres,  along  with  the  aerosol  models, 
were  used  to  determine  specific  amounts  of  a constituent  at 
a particular  altitude.  The  amounts  were  converted  Into 
scallxig  factors  for  a one-km  path  length  by  multiplying  each 
amount  by  one  km.  The  scaling  factors  were  used  for  the 
determination  of  cross  sections  from  transmission  charts. 
Numerical  values  for  cross  sections  were  computed  for  one-km 
Intervals  from  0 through  30  km.  The  aerosol  attenuation 
cross  sections  are  assumed  Identical  for  all  geographical 
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Table  I 

U.  S.  Standard  Atmosphere  (Hef  8:94) 


U.S.  STANDARD  ATMOSPHERE,  1962 


Ht. 

(km) 

Pressure 

(mb) 

Temp» 

(°K) 

Densi^ 

(g/m^) 

Water  Vapor 
(g.'m^) 

Ozone 

(g/m^) 

0 

1.  013EKI3 

288.1 

1.  225E+03 

5.  9E+00 

5.4E-05 

1 

8.  986E+02 

281.6 

1.  lllE+03 

4.  2E+00 

5. 4E-05 

2 

7.  950E+02 

275.1 

1.  007E+03 

2.  9E+00 

5.  4E-05 

3 

7.  312E+02 

268.7 

9.  093E+02 

1.  8E  >00 

5.  OE-05 

4 

6. 166E+02 

262.2 

8. 193E+02 

1.  lE+00 

4.6E-05 

5 

5. 405E-K)2 

255.7 

7.  364£-K)2 

6.  4E-01 

4. 5E>05 

6 

4.  722E+02 

249.2 

6.  601E+02 

3.  8E-01 

4.  5E-05 

7 

4.  lllE+02 

242.7 

5. 900E+02 

2.  lE-01 

4. 8E-05 

8 

3.  565E+02 

236.2 

5.  258E+02 

1.  2E-01 

5.2E-05 

9 

3.  080E+02 

229.7 

4.  671E+02 

4.  6E-02 

7.  lE-05 

10 

2.  650E+02 

223.2 

4. 135E+02 

1.  8E-02 

9.  OE-05 

11 

2.  270E+02 

216.8 

3.  648E-I-02 

6.  2E-03 

1.3E-04 

12 

1.  940E+02 

216.6 

3. 119E+02 

3.  7E-03 

1.6E-04 

13 

1.  658E+02 

216.6 

2.  666E+02 

1.8E-03 

1.  7E-04 

14 

1.  417E+02 

216.6 

2.  279E+02 

8.  4E-04 

1.9E-04 

15 

1.211E+02 

216.6 

1.  948E-!-02 

7.  2E-04 

2.  lE-04 

16 

1.  035E+02 

216.6 

1.  665E+02 

6.  lE-04 

2.  3E-04 

17 

3.  850E+01 

216.6 

1.  423E+02 

5.  2E-04 

2.  8E-04 

18 

7.  565E+01 

216.6 

1.216E+02 

4.  4E-04 

3.  2^-04 

19 

6.  467E+01 

216.6 

1.  040E+02 

4.  4E-04 

3.  5E-04 

20 

5.  529E+01 

216.6 

8.  891E+01 

4.  4E-04 

3.  8E-04 

21 

4.  729E+01 

217.6 

7.  572E+01 

4.  8E-04 

3.  8E-04 

22 

4.  047E+01 

218.6 

6.451E-I-01 

5.  2E-04 

3.  9E-04 

23 

3.  467E+01 

219.0 

5. 500E+01  . 

5.  7E-04 

3.SE-04 

24 

2.  972E+01 

220.6 

4.  694E+01 

6.  lE-04 

3.  6E-04 

25 

2.  549E+01 

221.6 

4. 008E+01 

6.  6E-04 

3.  4E-04 

30 

1. 197E+01 

226.5 

1.  841E+01 

3.  8E-04 

2,0E-04 

35 

5.  746E+00 

236.5 

8.  463E+00 

1.  6E-04 

l.lE-04 

40 

2.  871E+00 

250.4 

3.  996E+00 

6.  7E-05 

4.  9E-05 

45 

- 1.  491E+00 

264.2 

1.  966E+00 

3.  2E-05 

1.7E-05 

50 

7.  978fi-0l 

270.6 

1.  027E+00 

1.  3E-05 

4.  OE-06 

70 

5.  520E-02 

219.7 

8.  754E-02 

1.  5E-07 

8. 6E-0C 

100 

3.  008E-04 

210.0 

4.  989E-04 

1.  OE-09 

4.  3E-11 
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Table  II 

Tropical  Atmosphere  (Hef  8:3) 


TROPICAL 

Ht. 

(km) 

Pressure 

(mb) 

Temp. 

(°K) 

Density 

(g/m^) 

Water  Vapor 
(g/m3) 

Ozone 

(g/m3) 

0 

1.  013E+03 

300.0 

1.  167E+03 

1.  9E+01 

5.  6E-05 

1 

9.  040E+02 

294.0 

1.  064E+03 

1.  3E+01 

5.  6E-05 

2 

8.  050E+02 

288.0 

9.  689E4-02 

9.  3E+00 

5.  4E-05 

3 

7.  150E+02 

284.0 

8. 756E+02 

4. 7E+C0 

5.  lE-05 

4 

6.  330E+02 

277.  0 

7.  951E+02 

2.2E+00 

4.7E-05 

5 

5.  590E+02 

270.0 

7.  199E+02 

1.  5E+0C 

4.  5E-05 

6 

4.  920E+02 

264.0 

6.  SOlE-t-02 

8.  5E-01 

4.3E-05 

7 

4.  320E-f02 

257.0 

5.  855E+02 

4.7E-01 

4.  lE-05 

8 

3.780E+02 

250.0 

5.  258E-t-02 

2.5E>01 

3.  9E-05 

9 

3.  290E+02 

244.0 

4.  708E+02 

1.2E-01 

3.  9E<05 

10 

2.  860E4-02 

237.0 

4.202E+02 

5.0E-02 

3.  9E-05 

11 

2.470E+02 

230.0 

3.  740E+02 

1.7E-02 

4.  lE-05 

12 

2.  130E+02 

224.0 

3.  316E+02 

6.  OE-03 

4.3E-05 

13 

1.820E+02 

217.0 

2.  929E+02 

1.  8E-03 

4.  5E-05 

14 

1.  560E-I-02 

210.0 

2.  578E+02 

1.  OE-03 

4.  5E-05 

15 

1.  320E+02 

204.0 

2.260E-t-02 

7.  6E-04 

4.7E-05 

16 

1.  llOE-i-02 

197.0 

1.  972E+02 

6.4E-04 

4.7E-05 

17 

9.  370E+01 

19S.  0 

1.  676E+02 

5.  6E-04 

6.9E-05 

18 

7.  890E-I-01 

199.0 

1.  382E-t-02 

5.  OE-04 

9.0E-05 

19 

6.  660E+01 

203.0 

1.  145E-t-02 

4.  9E-04 

1.4E-04 

20 

5.650E+01 

207.0 

9.  515E-t-01 

4.  5E-04 

1.  9E-04 

21 

4.  800E-I-01 

211.0 

7.  938E+01 

5.  lE-04 

2.4E-04 

22 

4.  090E-I-01 

215.0 

6.  645E+01 

5.  IE -04 

2.8E-04 

23 

3.  500E+01 

2 17.  0 

5.  618E+0 

5 lE-04 

3.2E-04 

24 

3.  OOOE+01 

219.0 

4.763E+01 

6.  OE-04 

3.  4E-04 

25 

2.  570E+01 

22 1.  0 

4.  045E+01 

6.  7E-04 

3.4E-04 

30 

1. 220E+01 

232.0 

1.831E-H)1 

3.  6E-04 

2.4E-04 

35 

6.  OOOE+00 

243.0 

8.  600E-t-00 

1.  lE-04 

9.2E-05 

40 

3.  050E+00 

254.0 

4.  181E-t-0U 

4.  3E-05 

4.  lE-05 

45 

1.  590E+00 

265.0 

2.  097E+00 

1.  9E-05 

1.  3E-05 

50 

8.  540E-01 

270.0 

1.  lOlE-t-00 

0.  3E-06 

4.  3E-06 

70 

5.790E-02 

219.0 

9.210E-02 

1.  4E-07 

8.  6E-08 

100 

3.  OOOE-04 

210.0 

5.  OOOE-04 

1.  OE-09 

4.  3E-11 

13 


r^wy?. 


omt/PH/^'5-7 


■seasonal  modela. 

The  cross  sections  calculated  are  average  cross  sec- 
tions over  a given  band  width.  The  average  cross  sections 
are  weighted  by  the  Flancklan  distribution  of  photons  given 
by  the  6000^  average  temperature  source.  The  band  widths 
were  determined  by  sensitivity  studies  which  will  be  ex- 
plained In  Chapter  IV.  Average  cross  sections  for  aerosol 
and  molecule  absorption  and  scatter  are  summed  to  determine 
an  average  total  cross  section  (attenuation)  for  a given 
altitude. 


The  transmission  charts  used  to  calculate  the  aerosol 
and  molecule  cross  sections  are  all  similar  to  Figure  3. 

The  variation  of  transmittance  with  respect  to  wavelength 
can  be  determined  by  use  of  a vertical  transmittance  scale 
associated  with  a set  of  scaling  factors  for  the  absorber 
amount.  In  order  to  use  the  charts  efficiently ^ It  Is 
necessary  to  trace  the  transmittance  scale  and  associated 
factors  on  to  some  transparent  paper.  Thus  by  displacing 
the  transmittance  scale  horizontally  according  to  wavelength 
and  vertically  according  to  the  appropriate  equivalent  sea 
level  quantity,  one  can  construct  a transmittance  spectrum 
for  any  band  Interval.  The  transmittance  spectmn  can  be 
converted  to  an  average  cross  section  for  a band  by  weighting 
the  transmission  spectrum  with  a Flancklan  distribution 
representative  of  the  photon  source. 
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The  cross  sections  are  determined  from  the  transmit- 
tance, which  is 

T • ezp(-(iR)  (10) 

Solving  (10)  for  p yields 

(11) 

H 

where  R is  equal  to  one  km  and  is  the  cross  section  in 
units  km*^« 

The  equivalent  sea  level  absorber  amounts  for  a given 
atmosphere  or  aerosol  distribution  are  given  in  charts  like 
the  one  in  Figure  4.  The  a^^  (aerosol)  on  the  vertical 
scale  of  Figure  4 is  the  ratio  of  number  density  of  aerosols 
at  a particular  altitude  to  the  number  density  of  aerosols 
at  sea  level  for  a visibility  of  23  km  (Ref  8:34).  Using: 
Figure  4,  a^  at  an  altitude  of  zero  km  for  a visibility  of 
23  km  is  one.  The  aerosol  scaling  factor,  A (km),  for  a 
one-km  path  is  therefore  one  km.  If  the  tremsmittemce 
scale  of  Figure  3 Is  placed  on  the  chart  so  that  the 
scaling  factor,  one  km,  rests  on  the  datum  line  (as  shown 
in  Figure  3),  then  the  transmlttemoe  can  be  read  off  as  a 
function  of  wavelength  for  aerosol  absorption  or  extinction 
(absorption  plus  scatter).  The  chart  is  moved  horizontally 
until  the  transmittance  scale  is  positioned  directly  above 
the  required  wavelength  and  read  off  the  transmittance 
where  the  scale  crosses  the  curve.  For  example,  in  Flgiure 
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the  scale  Is  set  to  the  right  of  30  microns  and  the 
transmittance  Is  approximately  0.998.  If  the  transmittance 
Is  required  for  a different  aerosol  scaling  factor,  for 
example,  A » 0.01  km  for  a l»km  path  at  a 20-km  altitude, 
the  transmittance  scale  Is  displaced  vertically  until  the 
scaling  factor  0.01  Is  coincident  with  the  datum  line.  The 
scale  Is  then  moved  horizontally  to  the  require  wavelength 
and  read  off  the  transmittance.  The  one-km  transmittance 
can  be  converted  to  cross  sections,  per  kilometer,  and 
averaged  over  the  appropriate  baM  Interval. 

Using  charts  for  aerosol  aM  molecule  absorption  and 
scatter,  average  cross  sections  can  be  determined  for  every 
kilometer  from  0 through  30  km.  The  sum  of  the  Individual 
scatter  and  absorption  cross  sections  will  determine  the 
total  cross  section  daflned  In  Chapter  II.  These  cross 
sections  will  be  average  cross  sections  over  a given  band. 

The  average  total  cross  sections  determined  will  con> 
slder  aerosol  absorption  from  0.6  to  4.0  microns  and.  aerosol 
scatter  from  0.2  to  4.0  microns.  The  only  atmospheric  mole- 
cule considered  In  this  report  from  0.2  to  0.7  microns  Is 
ozone.  Therefore,  all  molecular  absorption  In  wavelengths 
less  than  0.7  microns  will  be  determined  from  ozone  amounts. 
Molecule  absorption  for  wavelength  greater  than  0.7  microns 
will  consider  water  vapo^°  and  uniformly  mixed  gases  (C02i 
N2O,  CO,  CHji^,  and  O2}  amounts. 

The  overall  accuracy  In  cross  sections  Is  limited  by 
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graphical  techniques.  KcClatchey  points  out  that  the 
accuracy  of  the  transmittance  which  this  technique  provides 
is  better  than  10  percent.  The  curves  presented  in  the 
transmittance  charts  tend  to  overestimate  the  transmittance 
for  very  long  paths  and  underestimate  the  transmittance  for 
very  short  paths  (Bef  8:37). 

The  average  total  cross  sections  for  given  band 
intervals  at  one-km  increments  from  0 through  30  km  are 
given  in  Appendix  B. 
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IV,  Sensitivity  of  Thermal  Pluence 
Prom  an  Air  Burst 

Before  ^creloping  ^ scnputer  code  to  calculate  the 
virgin  fluence  from  an  air  burst,  several  sensitivity 
studies  are  accomplished  to  determine  the  significance  of 
the  different  parameters  and  assumptions  considered.  The 
goal  of  the  sensitivity  studies  Is  to  determine  a reliable 
method  for  calculating  thermal  fluence.  The  choice  and 
number  of  energy  bands  used  to  specify  the  atmorpherlc 
attenuation  coefficients,  the  relative  Importance  of 
scatter,  the  significance  of  source  (fireball)  temperature, 
and  the  Importance  of  height  of  burst  will  be  considered. 

Number  and  Choice  of  Energy  Bands 

Prom  an  examination  of  the  range  of  thermal  photon 
energies,  the  grouping  of  the  photon  energies  Into  ultra- 
violet, visible  and  Infrared  bands  Is  the  most  natural  way 
to  group  the  energy  distribution.  The  Infrared  band  Is 
divided  Into  band  numbers  I (0.7  - 1»2  microns)  and  II 
(1.2  - 4,0  microns)  to  separate  the  "near"  and  "far" 

Infrared  respectively.  Average  total  cross  sections  for 
each  band  are  determined  by  methods  explained  in  Chapter  III. 
This  band  grouping  Is  defined  as  the  "broad"  band  model 
since  each  band  spans  a broad  wavelength  Interval. 

An  examination  of  HcClatchey's  transmittance  charts. 
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for  example  Figure  3.  Indicates  that  In  the  ultraviolet 
and  visible  bands  the  cross  sections  are  slowly  varying 
and  treatment  as  a single  broad  group  should  be  reasonably 
aocu3r&te,  ho^fever  In  the  infrared  region  large  discontinui- 
ties exist  In  the  cross  section  variations  with  photon 
wavelength  within  these  broad  bands.  These  atmospheric 
"windows*  (regions  with  no  attenuation)  for  water  vapor 
absorption  cross  sections  are  averaged  over  when  the  Infra- 
red  region  Is  divided  Into  bands.  A "fine*  band  Infrared 
model  Is  determined  by  dividing  the  Infrared  band,  0.4  - 
4.0  microns,  Into  l4  separate  bands.  The  divisions  between 
the  bands  are  determined  by  a visual  Inspection  of  the 
transmittance  chart  for  water  vcpor  absorption,  as  can  be 
seen  In  Figure  5*  Divisions  are  chosen  at  wavelengths 
where  large  discontinuities  are  obvious  from  a visual 
Inspection  of  the  transmittance  chart.  The  band  numbers  and 
wavelength  limits  (In  microns)  are  given  In  Table  5 In 
Appendix  C.  To  compare  the  "broad”  Euod  "fine"  band  models, 
a series  of  virgin  thermal  fluenco  calculations  were  made. 

The  virgin  thermal  fluence  Is  calculated  using 
equation  (6)  on  page  5 of  Chapter  II.  These  calculations 
are  based  on  three  considerations.  The  first  Is  a thermal 
yield  of  100  kllotons.  Secondly,  a 6000°K  blackbody  source 
Is  used  for  the  distribution  of  energies  by  band.  Finally, 
the  average  cross  sections  consider  only  the  absorption  due 


ONB/PH/75-7 

> 

to  water  vapor.  The  results  of  the  calculated  fluence 
(cal  per  cm  ) at  various  horizontal  ranges  are  given  In 
Table  III. 

The  conclusion  drawn  Is  that  the  "fine*  band  grouping 
gives  consistently  larger  values  for  the  fluence  calculated 
at  various  ranges  than  the  "broad"  band  grouping.  At  10-km 
range  the  difference  Is  nearly  a factor  of  8;  2.3  versus 

j 2 

0.3  cal  per  cm‘*.  Two  cal  per  cm  will  Ignite  shredded 
newspaper  (Ref  6:332). 

Importance  of  Scatter 

Many  thermal  calculations  ignor  scatter  of  the  thermal 
radiation  (Ref  2,  3 7)*  and  treat  only  the  virgin  (or 

unreacted)  fluence.  Some  scatter  does  occur.  The  question 
Is,  how  Important  Is  Its  contribution? 

The  relative  Importance  Is  measured  here  by  computing 
the  virgin  or  direct  fluence  with  two  different  attenuation 
cross  sections;  absorption  only  and  absorption  plus  scatter. 
The  true  answer  must  be  In  between  the  fluences  computed 
with  those  two  cross  sections.  Absorption  only  assumes 
scatter  does  not  occur  and  Is  the  upper  limit.  Absorption 
plus  scatter  assumes  every  scatter  as  removal  and  the  photon 
can  never  reach  the  target  point  even  after  multiple  scatter, 
and  therefore  It  Is  a lower  limit.  The  fluence  calculations 
were  made  for  a 200  klloton  burst  with  the  burst  altitude 
St  5 km  and  the  receiver  at  0 km.  Cross  sections  were 
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Table  III 

Comparison  of  "Pine"  £uad  "Broad"  Band 
Infrared  Fluence  Calculations 
(Table  entries  are  energy  fluence  In  cal/cm^) 


distance  between  burst  euid  receiver  (km) 
^and  numbers 
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determined  from  data  in  Reference  8,  The  results  of  these 
calculations  are  shown  in  Table  IV.  The  results  of  the 
dual  calculation  are  illustrated  in  Figure  6. 

The  conclusion  drawn  is  that  scatter  is  important  in 
computing  fluence  at  heights  below  five  km.  While 
determining  cross  sections  for  absorption  and  scatter,  it 
was  noted  that  aerosol  scatter  tends  to  be  the  dominant 
interaction  in  the  visible  region  for  altitudes  below  five 
km.  The  concentration  of  aerosols  in  the  lower  atmosphere 
is  reflected  in  the  contribution  of  scatter  in  the  com- 
putation of  thermal  fluence. 


Temperature  of  Source 

The  temperature  of  the  source  of  thermal  radiation 
can  be  modeled  approximately  by  a 6000°K  blackbody.  This  ‘ 
fireball  model  assumes  the  total  thermal  yield  is  radiated 
at  this  average  temperature.  Figure  1 on  page  3 of  Chapter 
II  illustrates  the  fact  that  the  temperature  of  the  second 
pulse  rises  rather  rapidly  and  cools  steadily.  During  the 
second  pulse,  the  fireball  is  radiating  thermal  energy  at 
different  temperatures.  The  temperatures  of  interest  range 
from  3500®K  to  7500®K.  John  P.  Cahill  has  modeled  a nuclear 
air  burst  by  time  weighting  the  temperatures  of  the  source 
to  account  for  the  different  temperatures  at  which  the 
source  radiates  thermal  energy  during  the  second  pulse. 
Cahill  found  the  source  radiated  at  an  effective  temperature 
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Table  IV 

PluATice  with  Absorption  plus  Scatter 
and  Absorption  only  Cross  Sections 


Distance  between 
burst  and  receiver 
(km) 

Pluence  with 
absorption  and« 
scatter  (cal/cm^) 

Pluence  with 
absorption^only 
(cal/cm^) 

5 

11.66 

15.28 

10 

2.21 

3.51 

15 

.77 

1.49 

20 

.35 

.81 

25 

.18 

.51 

30 

.10 

.35 

Pig,  6,  Pluence  with  Absorption  plus  Scatter  and 


Absorption  only  Cross  Sections  (Data  from 
Table  IV). 

t 
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of  7500®K  about  43. 5J^  of  the  time,  6500°K  for  12.5^,  5500®K 
for  19.8J^,  4500°K  for  l4.8j^  and  3500®K  for  9.4j^(Ref  1:17). 
Using  Cahill's  weighting  factors,  the  thermal  fluence  was 
calculated  for  a co«altltude  burst  and  receiver  of  10  km 
and  a thermal  yield  of  100  kllotons. 

The  total  virgin  fluence  calculated  at  the  receiver 
varied  less  than  five  percent  between  the  weighted  tempera- 
ture and  the  6000^  average  temperature  tmder  the  same 
conditions. 

The  conclusion  Is  that  the  computed  fluence  Is  rela- 
tively Insensitive  to  the  source  temperature  over  the  energy 
range  considered. 

Height  of  Burst 

The  thermal  fluence  from  a nuclear  air  burst  Is  affected 
by  the  height  of  burst.  The  thermal  yield  of  an  air  burst 
Is  a function  of  both  the  height  of  burst  and  total  weapon 
yield.  For  air  bursts  (less  than  30  km), -the  thermal  yield 
Is  released  during  the  second  thermal  pulse.  Empirical 
relationships  for  the  thermal  efficiency  and  thermal  yield 
have  been  developed  by  the  Analysis  Division,  Air  Force 
Weapons  Laboratory  (AFWL),  Klrtland  AFB,  New  Mexico.  A 
graph  of  the  thermal  efficiency  curve  developed  by  AFWL  Is 
Included  In  Reference  3«  The  thermal  efficiency,  t.  Is 
given  by 

T - exp(<-.358  - .009H  ♦ 7.14  x 10*V)  x 2.30)  x 

exp((-1.25  X 10" V ♦ 6.42  X 10"®H^)  x 2.30)  (12) 
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where  H is  the  height  of  burst  (km).  The  thermal  yield 
(kilotons),  Y,  for  an  air  burst  Is  given  by 

Y (13) 

where  W Is  total  weapor  yield  in  kllotons.  The  results  for 
a 200-klloton  total  weapon  yield  are  shown  in  Figure  ?• 

The  percent  of  photons  in  the  ultraviolet  region  in- 
creases as  the  burst  altitude  increases,  resulting  in  an 
apparent  rise  in  the  average  surface  temperature  of  the 
fireball  (Bef  7:8).  The  transmission  of  thermal  radiation 
is  affected  by  changes  in  the  composition  of  the  atmosphere 
with  changes  in  altitude.  All  of  these  factors,  thermal 
yield,  percent  of  photons  in  ultraviolet  region,  and  com- 
position of  the  atmosphere,  influence  the  calculation  of 
fluence  from  an  air  burst. 

The  transmission  of  thermal  radiation  is  affected  by 
both  the  height  of  burst  and  the  height  of  receiver.  In 
general,  the  higher  the  burst  and  receiver,  the  fewer  the 
molecules  and  aerosols  available  for  scatter  and  absorption. 
One  exception  is  ozone,  which  Increases  to  a mazimiim  concen- 
tration at  approximately  25  km.  The  effects  of  water  vapor 
are  primarily  confined  to  altitudes  less  than  10- km.  The 
aerosol  distribution  varies  with  the  surface  visibility 
only,  in  the  lower  five  km  of  the  atmosphere. 

The  thermal  radiation  code,  TRAP,  developed  for  APWL, 
uses  power  fractions  to  represent  the  distribution  of  thermal 
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Percent  Thermal  Yield 
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energy  over  different  bands  (Ref  7:8).  The  power  fractions 

are  based  on  empirical  data  and  distribute  the  thermal 

energy  approximately  like  blackbody  radiators.  The  percent 

of  thermal  radiation  in  the  ultraviolet  region  increases 

with  the  height  of  the  burst.  At  18  km  the  distribution  can 

be  approximated  by  a 700C^K  blackbody  and  at  32  km  by  an 

8000®K  blackbody.  Below  15  km,  the  6000°K  blackbody  is  a 

close  approximation.  The  power  fractions  were  compared  to 

blackbody  distributions  by  the  use  of  BAND  tables  for 

Plancklan  distributions  (Ref  5)*  Appendix  C illustrates 

the  Increase  in  ultraviolet  with  height  as  modeled  by  the 

changes  in  blackbody  temperature. 

The  general  decrease  in  absorption  and  scatter  with 

height  is  reflected  in  the  fluence  computations  in  Appendix 

D for  a 200-kiloton  burst  with  distance  between  burst  and 

receiver  of  10  km.  The  fluence  calculated  with  the  height 

2 

5f  burst  at  5 km  in  problem  one  was  2.21  cal  per  cm  , and 
fith  the  height  of  burst  at  20  km  in  problem  four,  thr 

O 

,‘luence  was  4.21  cal  per  cm  . Part  of  the  increase  in 
fluence  is  due  to  Increase  in  thei*mal  yield. 

The  important  conclusion  with  respect  to  height  of 
burst  is  that  burst  height  changes  can  be  modeled  by  changing 
source  temperature. 
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V,  Virgin  Thermal  Pluence  Code 

The  findings  of  the  sensitivity  studies  discussed  fji 
Chapter  IV  and  the  cross  sections  determined  by  methods 
outlined  in  Chapter  III  were  used  to  develop  a computer 
code  for  calculating  the  virgin  thermal  fluence  from  a 
nuclear  air  burst.  The  code  is  written  in  FOBTRAN  IV  com- 
puter language.  Comment  cards  are  used  throughout  the 
program  to  define  terms  and  identify  steps  in  the  calcula- 
tions. The  input  pareuneters  are  weapon  yield  in  kilotons, 
height  of  burst  in  kilometers,  height  of  receiver  in  kilo- 
meters, distance  between  burst  and  receiver  in  kilometers 
and  surface  visibility  of  23  or  5 km.  The  energy  distri- 
bution fractions  and  average  attenuation  cross  sections 
are  read  into  the  program  at  execution  time  to  simplify 
the  program  and  make  the  code  flexible  if  different  distri- 
butions and  cross  sections  are  considered. 

Thermal  Yield 

The  thermal  yield  of  an  air  burst  is  a function  of 
both  the  height  of  burst  and  total  weapon  yield  as  mentioned 
in  Chapter  IV.  Equations  (i2)  and  (13)»  presented  in 
Chapter  IV,  are  used  in  the  code  to  compute  thermal  yield. 

Attenuation  Cross  Sections 

Average  total  attenuation  cross  sections  were  computed 
for  each  of  the  l6  bands  at  one-kilometer  increments  from 
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^ ^ 1 4- ft  Ik  were  infrared,  radiation 

Rands  numbered  1 to  »*er 

0 to  30  km.  Banas  «uu. 

B.M  nu.^r  15  -«  Ultraviolet  ^ 

,6  .as  visitla  radiation.  Cross  seotions  ror  toth  t 
Standard  Atnosphere  and  the  Tropical  Atnosphsre  usr. 

.tnsd.  Kaon  atnosphere  has  the  loner  5-tn  cross  see  i^s 
o»puted  for  the  23  and  5-lc  surface  visihiUties. 

cross  sections  are  reed  into  the  progran  on  data  ca  . 

n»..d  on  whatever  atnospherlc  cross  sec- 
Oalculations  are  based  on  wnaue 

ti„«  are  read  into  the  progran.  The  cross 

detemined  by  methods  explained  in  Chapter  III^ 

listed  in  Appendix  B for  both  the  B.  S.  Standard 

Atuosphere  and  the  Tropical  Atnosphere.  All  the  cros 

Atuospnero  , ^ ^ data  found  in  Beferenoe  8. 

sections  were  detemined  based  on  d 

Wten  the  burst  and  receiver  are  at  different  altitudes  t 
nvemge  cross  sections  for  the  layer  of 

burst  and  the  height  of  receive 
the  l6  bands. 


m.trlbutl<ffi  of  SomSS  SSEBZ 

;r;ntributlon  .f  souro.  energy  ie  a function  of 

Che  blaohbody  temperature  used  for  the  fireball  n^el 

The  6000»K  blacKbody  was  used  for  all  burst  heights  15 
cr  less.  The  7000»K  blackbody  was  the  nodel  for  a urs 
.t  25  Kcr.  bursts  between  15  and  25  Bn.  linear  inter- 
,lla!ion  was  used  fron  the  SOCOOg  energy  distribution  to  the 
7000»K  energy  distribution.  For  burst  heights  greater 
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than  25  km,  linear  Interpolation  was  used  from  a 7000^ 
distribution  at  25  km  to  an  8000^  distribution  at  35  km. 

The  rise  In  average  fireball  temperatures  was  used  to  model 
the  observed  Increase  In  ultraviolet  radiation  with  height. 

The  band  divisions  are  given  In  Appendix  C with  the  relative 
distribution  of  energy  for  a 6000°K,  7000®K  and  8000°K 
Plancklan.  The  distribution  of  the  source  thermal  energy 
in  terms  of  ultraviolet,  visible  end  infrared  Is  part  of 
the  code  output. 

^Calculation  of  Pluenoe 

After  the  thermal  yield  has  been  calculated,  the 
energy  distribution  has  been  determined  and  average  cross 
sections  for  the  layer  between  the  burst  and  receiver 
heights  have  been  calculated,  the  virgin  fluence  contributed 
by  each  of  the  16  bands  Is  computed.  The  distance  between 
the  burst  and  receiver  Is  an  Input  parameter  used  In 
equation  (6)  on  page  5*  Additional  factors  have  been  In- 
cluded to  convert  kilometers  to  centimeters  and  kilotons  to 
calories.  In  order  to  give  the  fluence  In  units  of  calories 
per  square  centimeter.  The  ultraviolet,  visible  and  Infra- 
red virgin  fluence  computed  are  separate  outputs  along 
with  the  total  virgin  fluence. 

Options  are  Included  to  print  out  the  average  cross 
sections  computed  In  the  calculation  of  fluence  based  on 
the  layer  between  burst  and  receiver.  The  energy  distribution 
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oonputed  based  on  the  height  of  burst  is  cuiother  prlnt-out 
option.  The  fluence  in  each  of  the  16  bands  is  also  an 
optional  print-out. 

Appendix  D contains  four  sample  problems  illustrating 
the  use  of  the  code  in  calculating  virgin  fluence  from  an 
air  burst.  Problems  one  and  two  compare  the  use  of  the 
TJ.  S.  Standard  Atmosphere  cross  sections  with  the  Tropical 
Atmosphere  cross  sections  for  a surface  visibility  of  23  km. 
The  difference  in  the  calculated  total  virgin  fluence  is 
due  to  the  difference  in  the  infrared  fluence,  which  is 
caused  by  the  additional  water  vapor  in  the  Tropical  Atmos- 
phere. Comparing  problems  two  and  three,  the  difference  in 
the  total  virgin  fluence  is  due  to  the  difference  in  the 
visible  and  Infrared  fluences  caused  by  the  reduction  of 
surface  visibility  from  23  to  5 km.  This  illustrates  the 
effect  of  a haze  layer  on  both  the  visible  and  Infrared 
fluences.  The  haze  layer  reduced  the  total  virgin  fluence 
more  than  the  difference  caused 'by  an  increase  of  water 
vapor,  as  shown  in  the  comparison  of  problems  one  and  two. 
Problem  four  demonstrates  the  increase  in  total  virgin 
fluence  caused  by  Increased  thermal  yield  with  Increased 
height  of  burst  and  decreased  attenuation,  due  to  fewer 
aerosols  and  molecules  at  increased  height  of  burst  and 
receiver.  Notice  that  in  all  four  problems  the  ultraviolet 
fluence  of  contributions  is  less  than  one  percent  of  the 
total  virgin  fluence. 

Appendix  A is  the  program  listing  for  the  fluence  code. 
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VI.  Results  and  Recommendat Ions 


Results 

A sensitivity  study  was  made  on  four  variables  used 
In  the  calculation  of  nuclear  thermal  radiation.  The 
results  were: 

1.  Number  and  choice  of  energy  bands:  Visible  and 

ultraviolet  regions  could  be  treated  accurately 
by  one  band  each,  but  the  Infrared  region  needs 

14  energy  bands  to  account  for  atmospheric  "windows" . 

2.  Importance  of  scatter:  Scatter  Is  Important  In 

ccHuputlng  fluence,  particularly  at  heights  below 
five  km. 

3.  Significance  of  source  temperature:  Allowing  source 

emission  temperature  to  vary  from  7500°K  to  3500°K 
In  five  time  steps  made  only  a five  percent  change 
In  fluence  as  conpared  to  a 6000^K  constant  tempera, 
ture  source. 

4.  Importance  of  height  of  burst:  Increases  In  source 

altitude  were  modeled  by  Increasing  blackbody 
soiu’ce  temperature  to  8000®K  at  35  km. 

A code  was  developed  to  calculate  virgin  thermal  fluence 
from  a nuclear  air  burst.  The  source  of  the  thermal  radiation 
was  considered  a blackbody  radiator.  Cross  sections  were 
determined  by  the  use  of  transmitt Ion  data  applied  to  the 
U.  S.  Standard  Atmosphere  and  the  Tropical  Atmosphere 
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oompllad  by  APCBL.  All  scatter  and  absorption  Interactions 
were  considered  loss  mechanisms  along  with  attenuation  by 
spherical  divergence. 

The  energy  distribution  and  cross  section  sensitivities 
were  considered  while  developing  the  code.  A "fine”  band 
structure  was  found  to  be  more  representative  of  the  Infrared 
range  than  a "broad"  band  structure.  Increases  In  the 
ultraviolet  radiation  range  with  Increases  In  height  were 
modeled  by  Increasing  the  blackbody  temperature  with  height. 
By  changing  cross  sections,  the  effects  of  season,  geographi- 
cal location,  and  haze  layer  on  the  calculated  virgin  fluence 
were  demonstrated. 

The  code  was  written  to  be  easily  understood  and  used 
with  only  a basic  knowledge  of  FORTRAN  IV  language  and  a 
brief  Introduction  to  thermal  nuclear  effects.  The  program 
was  prepared  for  use  on  the  GDC  6600  computer  with  a Scope 
3.4  compiler  version.  Core  memory  required  on  the  GDC 
system  Is  approximately  40,000  octal  words  and  run  times 
are  In  the  order  of  a few  seconds. 

The  accuracy  of  the  code  calculations  for  virgin  thermal 
fluence  are  within  the  range  of  the  values  given  by  Glasstone. 
A series  of  four  problems  were  run  using  a 200-klloton  burst 
and  a distance  between  burst  and  receiver  of  10  km.  All  of 
the  total  virgin  fluences  calculated  for  surface  visibilities 
of  23  km  ranged  from  two  to  five  calories  per  cm  . Glasstone 
gives  a value  of  approximately  three  calories  per  cm  for 
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air  bursts  less  than  20  miles  with  the  distance  between 

burst  and  receiver  of  10  km  (Ref  6:333).  The  fluences 

calculated  by  the  code  for  bursts  less  than  five  km  were 

2 

all  less  than  three  cal  per  cm  at  a distance  of  10  km. 

The  lower  fluence  values  from  the  code  could  be  attributed 
to  variations  in  atmospheric  models. 

The  code  does  not  consider  reflection  from  the  earth's 
surface  or  cloud  layers.  The  only  haze  layer  considered 
Is  the  ground  based  haze  In  the  flve-km  surface  visibility 
model.  Attenuation  due  to  rain,  snow  and  clouds  was  not 
considered. 

Recommendat Ions 

Some  areas  that  could  be  further  Investigated  Include 
the  following  Items.  In  the  future,  different  cross  sections 
could  be  determined  using  other  atmospheric  models.  Since 
the  aerosol  model  used  has  the  same  distribution  of  aerosols 
above  five  km  ard  aerosol  cross  sections  were  the  largest 
part  of  the  atmospheric  attenuation  for  a number  of  cal- 
culations, other  aerosol  models  could  be  tried  to  Illustrate 
their  effects  on  thermal  fluence.  The  reaction  of  ultra, 
violet  photons  with  ozone  could  be  examined.  There  Is  a 
possibility  of  the  photo-decomposltlon  of  the  ozone  In  the 
atmosphere  by  the  initial  photons  from  a burst  and  the 
subsequent  photons  penetrating  to  the  receiver  because  of 
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the  depletion  of  the  absorbing  ozone.  The  photo- 
deoomposltlon  rate  and  the  rate  of  recombination  to  ozone 
aQuilibrlum  could  be  examined. 
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Appendix  B 

Cross  Sections  Listing 

The  average  total  attenuation  cross  sections  (per  km) 
for  the  IS  bands  in  the  code  are  listed  in  three  rows  for 
each  altitude  from  0 through  30  km.  The  first  two  rows  are 
the  l4  infrared  bands.  The  third  row  lists  band  15,  ultra, 
violet,  and  band  I6,  visible.  To  the  right  of  each  three, 
row  listing  is  the  altitude  of  the  cross  sections  in  km. 
Pages  47  through  49  list  the  U.  S.  Standard  Atmosphere  Cross 
Sections  for  a surface  visibility  of  23  km.  Page  50  lists 
the  haze  layer  cross  sections  for  the  U.  S.  Standard  Atmos, 
phere.  Pages  $1  through  53  list  the  Tropical  Atmosphere 
cross  sections  for  a surface  visibility  of  23  km.  Page  54 
lists  the  haze  layer  cross  sections  for  the  Tropical 
Atmosphere.  The  average  total  attenuation  cross  sections 
were  determined  by  methods  outlined  in  Chapter  III  with  all 
data  used  from  Reference  8. 
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Appendix  C 
Table  V 

Flaacklan  Distributions  of  Energy  by 
Bands  and  Blackbody  Temperature 
(Table  entries  are  percent  of  thermal  energy  i .i 


Band 

Numbers 


%5 

16 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

Total 


Band  Limits 
(microns) 


- 0.4 


0.7  - 0.75 
0.75-  0.85 
0.85-  0.95 
0.95-  1.05 
1.05-  lol5 
1.15-  1.25 
1.25-  1.50 

1.50-  1.70 
1.70-  2.00 
2.00-  2.50 

2.50-  3.00 
3.00-  3.50 

3.50-  3.80 
3.80-  4.00 

2 - ■ 


6000®K 

i%) 

7000°K 

i%) 

13.8 

22.7 

37.5 

39.2 

4.9 

4.4 

8.3 

6.9 

6.5 

4.5 

5.1 

4.8 

3.9 

3.1 

3.3 

2.6 

5.6 

4.0 

2.9 

2.0 

2.7 

1.9 

2.4 

1.7 

1.2 

0.8 

0.6 

0.5 

{ 0.3 

0.2 

0.1 

0,1 

99.1 

99.4 
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Appendix  D 
Sample  Problems 

VIRGIN  THEPMftL  FLUEtJCE  PRORLEH  1 

INPUT  PARAMETERS... 

WEAPON  YIELD  IS  230.  KILOTONS 

HEIGHT  OF  BURST  IS  5.  KM 

HEIGHT  OF  RECEIVER  IS  0.  KM 

DISTANCE  BETWEEN  BURST  AND  RECEIVER  IS  10.  KM 

PREVAILING  SURFACE  VISIBILITY  IS  23.  KM 

CROSS  SECTIONS  ARE  U.S.  STANDARD  ATMOSPHERE 

OUTPUT  PARAMETERS., . 

AT  THE  SOURCE... 

IR  thermal  YIELD  IS  28. S2  KILOTONS 
VISIBLE  thermal  YIELD  IS  22. A7  KILOTONS 
UV  THERMAL  YIELD  IS  8,?5  KILOTONS 
TOTAL  THERMAL  YIELD  IS  59.63  KILOTONS 

AT  THE  RECEIVER... 

IR  VIRGIN  FLUENCE  IS  1.22  C.AL/CM2 
VISIBLE  VIRGIN  FLUENCE  IS  , 9S  CAL/CM2 

UV  VIRGIN  FLUENCE  IS  .01  CAL/CM? 

TOTAL  VIPGIN  FLUt  NPr  IS  2, 21  CAL/CM2 
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VIRGIN  THERMAL  FLUENCE  PROBLEM  2 

INPUT  PARAMETERS... 

WEAPON  YIELO  IS  203.  KILOTONS 

HEIGHT  OF  RURST  IS  5.  KM 

HEIGHT  OF  RECEIVER  IS  0.  KM 

DISTANCE  BETWEEN  BURST  AMD  RECEIVER  IS  10.  KM 

PREVAILING  SURFACE  VISI^TLITY  IS  23.  KM 

CROSS  SECTIONS  ARE  TROPICAL  ATMOSPHERE 

OUTPUT  PARAMETERS... 

AT  THE  SOURCE... 

IR  THERMAL  YIELO  IS  28.52  KILOTONS 
VISIBLE  THERMAL  YIELO  IS  ZZ»U7  KILOTONS 
UV  THERMAL  YIELO  IS  8.25  KILOTONS 
TOTAL  thermal  YIELO  IS  59.88  KILOTONS 

AT  THE  RECEIVER... 

IR  VIRGIN  FLUENCE  IS  1,05  CAL/CM2 
VISIBLE  VIRGIN  FLUENCE  IS  .98  CAL/CM2 

UV  VIRGIN  FLUENCE  IS  .Cl  CAL/CM? 

TOTAL  VIRGIN  FLUcNCE  IS  2.03  CAL/CM2 
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VIRGIN  THtRMAL  FLUEKTE  PROOLEM  3 

INPUT  PARAMETERS... 

WEAPON  YIELD  IS  200.  KILOTCNS 
HEIGHT  OF  ^URST  IS  5.  t-'M 
HEIGHT  OF  RECEIVER  IS  0.  KM 

DISTANCE  3-ZTWEEN  BURST  A?  0 RECEIVER  IS  13.  KM 
PREVAILING  SURFACE  VISIBILITY  IS  5.  KM 

CROSS  SECTIONS  ARE  U.S.  STANDARD  ATMOSPHERE 

OUTPUT  PARAMETERS... 

AT  THE  SOURCE... 

IR  THERMAL  YIELD  IS  28.52  KILOTONS 
VISIBLE  THERMAL  YIFLO  IS  22.47  KILOTONS 
UV  THERMAL  YIELD  IS  8.25  KILOTCNS 
TOTAL  THERMAL  YIELD  IS  59.88  KILOTONS 

AT  THE  RECEIVER... 

IR  VIRGIN  FLUENCE  IS  .72  CAL/CM2 
VISIBLE  VIRGIN  FLUENCE  IS  .42  CAL/CM2 

UV  VIRGIN  FLUENCE  IS  .OC  CAL/CM2 

TOTAL  VIRGIN  FLUENCE  IS  1.13  CAL/CM2 
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VIRGIN  THERMAL  FLUENCE  PROBLEM  4 

INPUT  PARAMETERS... 

WEAPON  YIELD  IS  20D.  KILOTONS 

HEIGHT  OF  BURST  IS  20.  KM 

HEIGHT  OF  RECEIVER  IS  30.  KM 

DISTANCE  BETWEEN  BURST  AND  RECEIVER  IS  10.  KM 

PREVAILING  SURFACE  VISIBILITY  IS  5.  KM 

CROSS  SECTIONS  ARE  U. S.  STANDARD  ATMOSPHERE 

OUTPUT  PARAMETERS... 

AT  THE  SOURCE... 

IR  THERMAL  YIELD  IS  26.38  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  25.61  KILOTONS 
UV  THERMAL  YIELD  IS  12.17  KILOTONS 
TOTAL  THERMAL  YIELD  IS  66.72  KILOTONS 

AT  THE  RECEIVER,.. 

IR  VIRGIN  FLUENCE  IS  2.25  CAL/CM2 
VISIBLE  VIPGIN  FLUENCE  IS  1.95  CAL/CM2 
UV  VIRGIN  FLUENCE  IS  .CC  CAL/CM2 

TOTAL  VIRGIN  FLUENCE  IS  4.21  CAL/CM2 
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